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APPROXIMATE ANALYSIS OF EFFECTS OF LARGE DEFLECTIONS AND INITIAL TWIST ON 
TORSIONAL STIFFNESS OF A CANTILEVER PLATE SUBJECTED TO THERMAL STRESSES 1 

By Richard R. Heldenfels and Louis R. Vostebn 


SUMMARY 

An approximate analysis of the nonlinear effects of initial 
twist and large deflections on the torsional stiffness of a cantilever 
plate subjected to a nonuniform temperature distribution is 
presented . The Von Kdrmdn large-deflection> equations are 
satisfied through the use of a variational principle . The 
results show that initial twist and applied moments can have 
significant effects on the changes in stiffness produced by non- 
uniform heating , particularly in the region of the buckling 
temperature difference , Results calculated- by this approximate 
analysis are in satisfactory agreement with measured- torsional 
deformations and changes in natural frequency, 

INTRODUCTION 

One of the structural problems of high-speed flight is the 
reduction of effective stiffness of structures due to the thermal 
stresses produced by aerodynamic heating. A reduction in 
torsional stiffness can be an important factor in aeroelastic 
problems as indicated in references 1 and 2. A similar 
reduction in stiffness produced by thermal stresses pre- 
sumably caused the flutter and failures of some structural 
models described in reference 3. A simple method for 
calculating the reduction in torsional stiffness of thin wings 
is presented in reference 4. In reference 5 the results cal- 
culated from a small-deflection plate theory are compared 
with experimentally determined changes in the torsional 
stiffness of a cantilever plate rapidly heated along the 
longitudinal edges. The theory used in reference 5 predicted 
the general effect of thermal stresses on the torsional stiffness, 
as indicated by measurements of torsional deformation and 
changes in natural frequency of vibration, hut overestimated 
the magnitude of the changes. 

The purpose of this paper is to present the results of an 
approximate analysis to show that the differences between 
theory and experiment noted in reference 5 are due to the 
nonlinear effects of large deflections and initial deformations 
not included in the small-deflection analysis. The analytical 
approach used to account for large deflections and initial 
deformations is presented, the general significance of the 
results is discussed, and calculated values are compared with 
the experimental data of reference 5, 
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1 Supersedes NACA Technical Note 4067 by Richard R. Heldenfels and Louie F. Vos teen, 1957. 


SYMBOLS 

plate length in ^-direction 
coefficients of series expansion for plate 
deflection 

half-plate width in ^-direction 
coefficient of stress function 

Ei ? 

plate flexural stiffness, -r^u\ £\ 

/X ) 

modulus of elasticity 

stress function defining stress distribution in 
plate 

selected functions of x and y 
values obtained from definite integrals 
torsional stiffness 

torsional stiffness of flat, unstressed plate 
incremental torsional stiffness 
initial incremental torsional stiffness before 
heating 

minimum incremental torsional stiffness 
nondimensional moment 
applied moment (positive in direction of 
positive twist) 
pressure 
plate thickness 
temperature 

particular temperature difference 
critical value of A T 
total plate deflection 
initial plate deflection 
coordinate axes 

coefficient of thermal expansion 
small dynamic perturbation 
exponential parameter of temperature dis- 
tribution function 
twist at plate tip 
initial twist at plate tip 
temperature ratio, ATjAT^ 
value of X corresponding to A(OJ) min 
Poisson's ratio 
density 


593 



594 


REPORT 136 1 — NATION AD ADVISORY COMMITTEE FOR AERONAUTICS 


cr x ,< 7 y normal stresses in plane of plate in x- 

and ^/-directions, respectively, positive 
for tension 

shear stress in plane of plate 
r time 

<p nondimensional twist 

(p { initial nondimensional twist 

a) circular frequency 

co 0 circular frequency of unheated perfect plate 

co j initial circular frequency 

co miH minimum circular frequency 

V 2 differential operator, 

b 4 d 4 b 4 

7 41 differential operator, ^+2 ^ dy i + ^d 

Double dot indicates second derivative with respect to 
time. 

ANALYSIS 

STATEMENT OF PROBLEM 

The studies presented herein are primarily concerned with 
the twist of a uniformly thick rectangular cantilever plate 
shown in figure 1. The equations are derived to consider the 
effects of initial twist, applied moments, and thermal stresses 
on the torsional deformations , and natural frequencies of 
plates. The analysis involves an approximate solution of 
the Von K&rrmin large-deflection equations which have been 
modified to include the effects of initial imperfections and 
nonuniform temperature distributions. The modified equa- 
tions are (from ref. 6) 


V 4 F=— Ea^T+E\ 


b a iob*u? 
bx 2 by 2 


w / bhVjV b*u? f b^r 
t \bx by/ bx 2 by 2 _ 


p_. _t_ / b 2 F b*w 
"ITD \by 2 bx 2 ' 


b»Fb^ b'F bho > 
bx 2 by 2 bxbybxbyy 


where F is the stress function such that 


Cx ~W 



Equation (la) is the differential equation for compatibility 
in the plane of the plate and equation (lb) is the differential 
equation for equilibrium of forces acting perpendicular to the 
plane of the plate. The solution of equations (1), subject to 
the proper boundary conditions on the deflections and 
stresses, describes the behavior of the plate under the applied 
loads p and the temperature distribution T. For the purpose 
of this paper the assumption is made that no external loads 
are applied in the plane of the plate so that the inplane 
stresses defined by the stress function F are zero at the free 
boundaries (stress-free edges). Equations (1) apply to both 
the dynamic and the static problem if p includes the dynamic 
loads as well as static loads. 

APPROXIMATE SOLUTION OF EQUATIONS 

Exact solutions of the Von K&rmdn large-deflection equa- 
tions are difficult to obtain, but approximate solutions can be 
obtained to any desired accuracy by several methods. For 
the analysis herein equations (1) will be satisfied through the 
use of the following variational principle: 

•Off { _ 


SfLW) -(w) j -2 ^ (S7 aT~5T Tf)/ 

JJ tocT (^ + w) dxdy ~§§ p ( y > —' w *) (Ix d y)= 0 (3) 

This principle is essentially a particularization of one given 
by Reissner (ref. 7) and also has been modified to include the 
effects of initial deformations and a nonuniform temperature 
distribution. The deflection w must satisfy the geometrical 
boundary conditions, that is, the conditions imposed on slopes 
and deflections. Variation of equation (3) with respect to F 
yields equation (la) and the associated natural boundary 
conditions, and variation with respect to w yields equation 
(lb) and the associated natural boundary conditions. 

The following assumptions for stress, deflection, and tem- 
perature are made : 

F=Of l ' 

w=bef» 

, r ( 4 ) 

«><=* >9 if 3 


Figure 1. — Dimensions and coordinate system of cantilever plate. 
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where the stress coefficient C and the tip twist 8 are to be 
determined by means of the vaiia tional principle* The func- 
tions/^/^ and/ 3 , the initial tip twist 6 if and the temperature 
coefficient AT are presumed to be known. Selection of these 
quantities is discussed in appendix A. It should be noted 
here that, although the temperature distribution is considered 
to vary with time, the shapes of the plate deflection and of 
the stress function are assumed to remain constant during 
heating or loading. 

The unknown coefficient of the stress function C is obtained 
by substituting equations (4) into equation (3) and taking the 
variation with respect to C . The result is 


and 8 is the Dirac delta function. 

If the value of C from equation (5) is substituted into 
equation (9) and the terms are rearranged, the resulting 
expression is 
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which gives the relationship between twist, applied moment, 
inertia loading, and temperature. 


C=—ExAT (~^+Eb t (8 1 —6i T ) (£) 


(5) 


where 

h 




m 


ffira (M, o*/, t'Mj,\‘_„ _a%_ Ml Mil 

5 JJ 2L&!/ 1 \bxj \by J bxby bx d?/J 1 ' 

The relation between the twist, the load, and the tempera- 
ture is obtained by performing the variation on equation (3) 
with respect to the twist 6 with C held constant. The 
resulting relationship is 


Db 2 (6— dijl+^tCPdls— bjjpfi dx 0 


(9) 


where 

4 


< io > 


The pressure loading p now will be considered to consist of 
a static loading p t and a d 3 mamic or inertia loading pt w so that 


■p=zp t — pt w 


(ii)' 


For this analysis the lateral static load p t will be restricted 
to two equal concentrated loads P applied at the corners 
x=a 7 y=±b in such a way as to form a couple about the 
z-oxis. Then 

JJVj dxdy=^I s —ptb0l e (12) 

where 

8(r—a)[8(y—b)—8(y+b)]f 2 dx dy (13) 

/e=JJ/s 2 <fa<fy (14) 

M=2bP (15) 


DEFINITION OF PARAMETERS 

It is convenient to define certain parameters which can be 
obtained from equation (16) when the large-deflection effects 
are neglected and the initial twist 8 t is 0. Then, from 
equation (16) 
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^4=~4-P^*4d 
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(17) 


If the plate is not. vibrating and no loads are applied, the 
critical, or buckling, temperature difference is obtained from 
equation (17) as 


, D PE 

2 Eat LJ Z 


(18) 


Similarly, when the plate is not vibrating and the tempera- 
ture is uniform ( AT=0 ) but the load is acting, the moment- 
twist relation is given by 


Ma 

\QJ) 0 


(19) 


where ( OJ) 0 is the torsional stiffness of the flat plate given by 

( 20 ) 


{GJ) 0 =2abwfy 

-LS 


In the absence of heating and loading the frequency may be 
found to be 


* Z>4 

w= 7tZ 


( 21 ) 


When equations (18), (19), (20), and (21) are substituted 
into equation (1 6) the result may be written as 

e -°'=Mkrb i+ ^J- H *- en ^W. (22) 

If a “nondimensional” twist is defined so that 


<P=6^ : 


2Etb 2 (4) 2 
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and a temperature ratio 


X— 


AT 
A Tn 


(25) 


equation (22) may be written simply as 

<p—<Pi=m+\<p— j <p 

C^o 


(26) 


APPLICATION TO STATIC PROBLEM 

The relationship between twist, moment, and temperature 
is given by equation (26) when the dynamic term is 0. 

COo 

Equation (26) then becomes 


cp— (pi—m+Xip— <p * 0 (27) 

The incremental stiffness of the plate, defined as the rate of 
change of moment with respect to twist, is given by 

A(OJ) a MI 

(6J)o i0J) o 

or, from equation (27), 

A ( QJ ) 1 *\ 1 0 9 2 /OQ\ 

L-X+V-w ( 28 ) 

APPLICATION TO DYNAMIC PROBLEM 

In order to determine the effect of temperature, moment, 
and initial twist on the natural frequency of torsional 
vibration, the quantity (p in equation (26) is replaced by 
<P~\~ € sin o>t where <p is considered to be the static solution 
obtained from equation (27) and e represents a small dy- 
namic perturbation about the static equilibrium position. 
Subtracting equation (27) from the perturbated relation, 
neglecting higher order terms in e, and dividing by e sin cor 
yields 

i =x-v+ *’ + © 

The frequency may then be written in terms of the twist 
and temperature as 


Q^l-X+3^-^ • (29) 

which is identical to the incremental stiffness (eq. (28)). 
This identity results from assuming similar deflection mode 
shapes for thermal buckling, twist due to applied moment, 
and torsional vibration and would not be expected to apply 
if the mode shapes involved were different. 

RESULTS AND DISCUSSION 

GENERAL RESULTS OF EQUATIONS 

The behavior of the plate as determined by equation (26) 
is discussed in the following sections for several combinations 
of conditions. In most cases the calculations include large 
values of initial twist and applied moment and have been 
carried well into the region where the temperature difference - 
AT exceeds the critical value. Although some of the results 


may be beyond the range for which the analysis is accurate, 
these results have been presented to illustrate the trends 
indicated by the equations despite the fact that they may 
not be quantitatively correct in some regions. 

Twisting duaio an applied moment. — The twisting of the 
plate due to an applied moment is given by equation (27) if 
the temperature ratio X is set equal to 0. The results from 
this expression are presented in figure 2 for various values of 
the initial twist. 

In figure 2, small-deflection results would plot as lines at 
45° to the coordinate axes. The large-deflection results 
become increasingly different as the moment or initial twist 
is increased and indicate that the plate becomes substan- 
tially stiff er as the twist is increased. 

Buckling due to nonuniform heating, — The twist of the 
plate (assuming that the buckling mode is a twisting action) 
is given by equation (27) when the moment m is 0; for the 
initially flat plate <pi would also be 0. The initially flat, or 
perfect, plate begins to deform only after the critical tem- 
perature difference is reached, whereas, as indicated by 
equation (27), the plate with initial twist begins to deform 
immediately upon heating. The results obtained from 
evaluation of equation (27) are given in figure 3 for several 
values of the initial twist. 

If the initial twist is large, a plot of twist against temper- 
ature ratio, like figure 3, does not give an accurate indication 
of the buckling temperature. Only if the initial twist is 
small is there a definite knee in the curve as the buckling 
temperature is approached, but this knee occurs below tho 
buckling temperature of the perfect plate. 

Combined action of applied moment and nonuniform 
heating. — Equation (27) applies directly to this case but 
may be more conveniently written as 

*-*/--XV+^-V 2 )==0 (30) 



Figure 2. — Calculated twist as a function of applied moment for 
' several values of the initial twist. 
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Figure 3. — Calculated twist as a function of temperature ratio for 
several values of initial twist. 


where 

(31a) 

and 

X'=\— (m 2 +2m<Pi) (31b) 

Equation (30) now has the same form as equation (27) when 
m is set equal to 0. The results plotted in figure 3 then apply 
to this case also if <p t is replaced by <p f and X by X'. The 
moment then effectively acts the same as an initial twist 
if the temperature ratio is reduced by the quantity m 2 +2m<pi. 

In figure 4 the relation between moment and twist has 
been indicated for two values of initial twist and several 
values of X, These curves show the characteristic increase 
in stiffness as the twist increases but also show a reduction 
in stiffness as the temperature ratio increases. These changes 
in stiffness are examined further in the following section. 
The portions of the curves where the slope is negative have 
been shown as dashed lines and are regions of unstable 
equilibrium which would not exist in the physical problem. 
The unstable portion of the curve exists whenever 
3?> 2 >L 

Frequency and incremental stiffness, — As has been noted, 
the square of the frequency ratio (eq. (29)) varies in the 
same manner as the incremental stiffness (eq. (28)) and, 
therefore, any of the following discussion pertaining to 
stiffness applies directly to the square of the frequency. 
Also, the figures which are presented for the incremental 
stiffness have been labeled with the square of the frequency 
ratio as well as with the stiffness ratio. 

The variation in incremental stiffness with temperature 
given by equation (28) is shown in figure 5 for several values 
of initial twist when m is equal to 0. The results show that, 
when the perfect plate is heated, the stiffness decreases 
linearly as the thermal stresses develop and becomes zero 
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Figure 4. — Calculated twist as a function of applied moment for 
several combinations of temperature ratio and initial twist. 



Figure 5. — Calculated incremental stiffness as a function of tempera- 
ture ratio for several values of initial twist. 

when the buckling temperature difference is reached 
Further heating causes the plate to twist and the twisting 
in turn leads to an increase in stiffness at a rate twice that 
of the initial decrease. 

If the plate has an initial twist, the heating causes the 
twist to increase as the heating progresses. This twist leads 
to an increase in stiffness that tends to counteract the reduc- 
tion produced by thermal stresses and, as a result, the 
stiffness first decreases and then increases without going 
to zero. The points of minimum stiffness occur at temper- 
atures lower than the buckling temperature. The locus 
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of points of minimum stiffness has been indicated in figure 
5 by the dashed fine. If the initial twist is sufficiently 
large, the stiffness does not decrease but increases with 
heating. 

When a moment is applied to the plate, figure 5 will 
indicate the incremental stiffness of the plate if <pt is replaced 
by <p' and X by X' (as defined by eqs. (31a) and (31b)). 

The locus of points where the incremental stiffness is a 
minimum is given by the equations 


MOJ) ntn '( Vi +m\™ 

~mr~\ 4 ) 

(32) 

Xmf*=l <Pt Oi ^ ) 

(33) 

«?J)o 2(1 Vi) 

(34) 


Equations (33) and (34) show that the nonuniform heating 
will increase the stiffness of the plate if the initial twist is 
geater than that given by 

The initial incremental stiffness of the plate is given by 
equation (28) when X is 0. The ratio of the minimum to 
the initial incremental stiffness therefore may be written as 



Initial incremental stiffness as a funotion of initial twist for several 
values of applied moment. 


A (QJ)min 
A ((?J) i 


1+3 


(35) 


Some results from equation (35) have been plotted in figure 
6. This figure indicates that small changes in the initial 
twist may cause large changes in the minimum incremental 
stiffness. The effects of applied moment are again similar 
to those of initial twist and certain combinations of the two 
can lead to drastic stiffness changes. For negative values 
of <p t the curves would be similar, with the stiffness ratio 
going to zero whenever pi— — w?. 

COMPARISON WITH EXPERIMENT 

Results calculated with the previously derived equations 
are compared in the following sections with the experimental 
results reported in reference 5. The results given in reference 
5 are for a square cantilever plate which was heated along 
the two longitudinal edges by carbon-rod radiators. Typical 
temperature histories of points on the heated edge and on 
the longitudinal center line are given in figure 7. Heat 
was supplied to the plate edges for 16 seconds; then the 
plate was allowed to cool. The deformations of the plate 
under the influence of nonuniform heating were determined 
for the conditions of no load and applied positive and 
negative tip moments. The changes in natural frequency 
of the first torsion mode were also measured during a heating 
test. J 

The expressions used in reference 5 for stresses, deflec- 
tions, and temperature are retained except that herein only 



Figure 7. — Measured temperature histories for a heated edge and 
the center line of the plate test specimen. 


three terms are used for the deflections. In nondimensional 
form these expressions are 


and 




(36) 

(37) 

(38) 
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Tiic function fi satisfies the condition of zero stress on the 
free edges and / a specifies zero slope and deflection along 
the root. The two undetermined coefficients of equation 
(37) are established from the small-deflection buckling 
analysis and carried through the large-deflection analysis as 
constants. Their values are given in appendix A along with 
data on the influence of the number of terms in the deflec- 
tion function on the accuracy with which the buckling 
temperature difference, natural frequency, and ratio of twist 
to an applied moment can be calculated. The exponential 
parameter f of the temperature function is adjusted to 
approximate the variation in the temperature distribution 
during a test. 

The numerical evaluation of the integrals and related 
functions required for comparison of theory and experiment 
are given in appendix B. 

Initial plate shape. — The initial shape of the plate of 
reference 5 was measured and is indicated in figure 8, In 
the selection of a value of initial twist, the higher order 
shapes (which are unlikely to have much influence on the 
twisting) have been ignored. The free comers were con- 
nected by a straight line in order to obtain a value of initial 
twist of 0<=O.35° which was used for comparison of theory 
and experiment. 

Twist due to an applied moment. — The deformations 
resulting from heating for three values of applied tip moment 
are presented in figure 9 and compared with curves calculated 
by use of a value of 0 f =O.35°. The agreement between 
theory and experiment is satisfactory, although the theory 
overestimates the twist in the vicinity of the maximum 
temperature difference (about 16 seconds). No theoretical 
results are presented past 20 seconds because, beyond this 
time, the actual temperature distribution cannot be repre- 
sented very well by the one-parameter temperature function. 

Another comparison of measured and calculated deforma- 
tions is shown in figure 10 where the abscissa is the tempera- 
ture difference instead of time. Because the shape of the 
temperature distribution changes with time, the portion of 
the curve for decreasing AT does not retrace the heating 
portion. 

Natural frequency. — The changes in natural frequency 
during a heating test are shown in figure 11 and compared 



y/in. 


Pi o dub 8. — Initial shape of plate test specimen. 
526097 — 60 39 



Figuee 9. — Comparison of measured and calculated values of twist 
as a function of time -for the rapidly heated plate for three values of 
applied moment. 



Fig dee 10. — Comparison of measured and calculated values of twist 
as a function of temperature difference for three values of the ap- 
plied moment. 
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FiGtrRE 11. — Comparison of measured and calculated values of fre- 
quency ratio as a function of time for the rapidly heated plate. 


with a calculated curve for 0,= 0.35° and with the small- 
deflection results of reference 5. The results calculated by 
the use of the large-deflection equations are in good agree- 
ment with experimental values and account for the frequency 
increase measured in the vicinity of the buckling temperature. 
A substantial improvement is noted over the predictions of 
small-deflection theory. In figure 12 the frequency ratio 
is shown as a function of temperature difference. 



Figure 12. — Comparison of measured and calculated values of the 
frequency ratio as a function of temperature difference. 


CONCLUDING REMARKS 

An approximate analysis of the effects of initial twist and 
large deflections on the torsional stiffness of a cantilever 
plate subjected to nonuniform heating shows that for a 
perfectly flat plate the effective stiffness, and thus the 
torsional frequency, decreases with increasing thermal 
stress, just as predicted by small-deflection theory, and 
goes to zero when the buckling temperature difference is 
reached. Beyond the buckling temperature, however, the 
stiffness, and thus the frequency, increases as the plate 
twists. If the plate has an initial twist, it begins to deform 
immediately upon heating and the stiffness decreases in 
much the same way as that of the perfect plate. The incre- 
mental stiffness of the initially twisted plate, however, 
reaches a minimum greater than zero before the theoretical 
buckling temperature difference is reached; further heating 
then increases the stiffness. The minimum incremental 
stiffness is a function of the initial twist and. if the initial 


twist is sufficiently large, no reductions in stiffness are 
obtained and nonuniform heating then always increases 
the stiffness. 

These results for stiffness changes associated with the 
torsional frequency are also applicable to the stiffness of 
the plate with respect to small changes in the applied 
moment. In this case the applied moment has an effect 
similar to the effect of an initial twist. If the applied moment 
exactly counteracts the initial twist, the plate behaves in 
much the same way as the perfect plate. 

Calculated results were compared with available experi- 
mental data and were found to be in satisfactory agreement 
in view of the approximate nature of the calculations. 

/ 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 15, 1957. 





APPENDIX A 

DISCUSSION OF MODAL FUNCTIONS USED IN THE ANALYSIS 


la the analysis of the torsional deformations of the 
cantilever plate by small-deflection theory, a deflection 
function containing six terms antisymmetries! in the co- 
ordinate y was used. The function (eq. (8) in ref. 5) is 

w^A^+Au^+Ai^+Ai^+A^+Au^y 1 (Al) 

The effect of various combinations of terms on the cal- 
culated critical temperature difference, torsional frequency, 
and moment-twist relation was investigated for a square 
cantilever plate (a/b= 2) by use of the small-deflection theory 
of reference 5. The following table shows the results obtained 
by starting with a single term A n and then progressing 
through other combinations to six terms: 


Terms used 

A TV 
Ki 

lCx 

M 

e 

-A\j - . 

103.9 

67.7 

60.4 

57.2 

60.0 

127.0 

77.5 

71.3 

77.4 
70.1 

0.324X10’ 

.278 

.273 

.278 

.273 

Atf f Aj} __ _ . 

A J2, An* 

Aft, At2i An, T - 

Ait, An* Ait* Au, Ati , Am 


where Ki and K 2 represent constants that include the plate 
dimensions and material properties. 

Inasmuch as the tabulated values were obtained from an 
application of the Rayleigh-Ritz procedure, the lowest 
value is the most accurate. The use of six terms improves 
the accuracy less than 2 percent over three terms and, 
consequently, is not worth the extra complication. The 
large-deflection analysis thus can be made by use of only 
three terms. The deflection function is nondimensionalized 
for this analysis and expressed in the form 


h — 


(P'WSI 


fjj 

l-}“ A? 


(A2) 



Figueb 13. — Comparison of deflected shapes for thermal buckling, 
torsional vibration, and twist due to an applied moment as given by 
sm all-deflection theory. 


data show that the deflection shapes are all approximately 
the same. Then, so long as the calculations are used to 
indicate changes from the initial conditions, any of the 
modes indicated should be satisfactory. 

In addition to the deflection function of equation (A2), 
a stress function aUd a temperature function are required in 
the analysis and have been selected to correspond to those of 
reference 5. For convenience, the functions are expressed 
in nondimensional form as follows: 


The relative values of these coefficients for the conditions 
given in the preceding table were obtained from the small- 
deflection analysis of reference 5 for the square cantilever 
plate (a/b= 2) with the following results: 



(A3) 

(A4) 


Condition 

1 

Ai 

As 

1+At+Ai 

1+At+At 

1+At+As 

Thermal bucklinp (\=1, 

■4,81 

-6.75 

2.14 

Torsional vibration (X=0, £c/t*=l)_ 

3.63 

— 4. 06 

1.42 

Applied moment 

2.38 

-2.01 

0.66 

Applied moment (X^0.6) 

3.10 

-3.20 

1.10 


The variation of / 3 with x/a for these four conditions has 
been plotted in figure 13. The table indicates a wide varia- 
tion in the relative values of the coefficients but the plotted 


The exponent f must be selected to describe best the meas- 
ured temperature distribution at the time of interest. Note 
that equation (A4) requires that the specified temperature 
difference AT be the difference between temperatures at the 
heated edges and the longitudinal center line. The func- 
tion/! has been chosen so that the stresses vanish on the free 
boundaries. Along the root (r=0) the function requires 
that the plate be free to expand parallel to the y-axis. 
Although this requirement does not indicate a “built-in ,> 
condition, it probably resembles the actual test condition. 
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APPENDIX B 


EVALUATION OF INTEGRALS AND RELATED FUNCTIONS 


The integrals indicated in the analysis have been evaluated 
in terms of the plate dimensions and the nondimension al 
functions given in appendix A. The integrals are as follows : 

'-ri 

f)'IKi)‘+*(lH 

(Bl) 

* - f *■ t ( W) W ) <“» 

128 (4 A t * 16 AS At 16 A 2 AtA s \ 
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For the square cantilever plate (a/b= 2) and the values of 
Ai and A 2 given in appendix A for the thermal-buckling 
mode, the integrals I Y to I 4 reduce to the following when 
Poisson's ratio /x equals 0.33: 
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The integrals / c and h have not been included but are 
discussed subsequently. 

The critical temperature difference given by equation (18) 
then becomes 
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The twist and applied moment are related to their cor- 
responding nondimensional quantities by equations (23) and 
(24) and may now be expressed as 

<P= 0.035 (B12) 
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Use of the thermal-buckling mode yields incorrect values 
for the natural frequency and the twist due to applied 
moment. In addition, the measured values differ slightly 
from the correct calculated values as a result of the im- 
perfect clamping at the root of the cantilever plate. For 
these reasons, the measured initial values of torsional fre- 
quency and twist due to an applied moment are used in the 
calculations and the changes produced by the heating cycle 
are calculated from the following equations : 


(Pi+m , 2 2 

A =1 V<P — <Pi 

<P 


(B14) 

(B15) 


The integrals I h and 7 0 are needed only to calculate the 
initial natural frequency and twist due to applied moment 
and, thus, have not been evaluated for the mode shape of 
thermal buckling. 

In the calculations, results of which are presented in 
figures 9, 10, 11, and 12, the following quantities were used: 

#= 10.6X10® psi 
a= 12. 8X10“® °F“ 1 
^=0.33 
£=0.25 in. 
a=20 in. 

6=10 in. 

0,=O.35° 
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